Crossover from the Josephson effect 
to bulk superconducting flow 



Fernando Sols and Jaime Ferrer 
Departamento de Fisica de la Materia Condensada, C-XII 
Universidad Autonoma de Madrid, E-2804-9 Madrid, Spain 

Abstract 

The crossover between ideal Josephson behavior and uniform superconduct- 
ing flow is studied by solving exactly the Ginzburg-Landau equation for a 
one-dimensional superconductor in the presence of an effective delta function 
potential of arbitrary strength. As the effective scattering is turned off, the 
pairs of Josephson solutions with equal current evolve into a uniform and a 
solitonic solution with nonzero phase offset. It is also argued that a micro- 
scopic description of the crossover must satisfy the self-consistency condition, 
which is shown to guarantee current conservation. The adiabatic response to 
an external bias is briefly described. The ac Josephson effect is shown to break 
down when the external voltage is applied at points which are sufficiently far 
from the junction. 
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I. INTRODUCTION 



The Josephson effect between two weakly coupled superconductors and the steady flow 
of supercurrent in a perfect lead constitute the two main paradigms of superconducting 
transport. Both regimes can be viewed as the limits of a general scenario in which Cooper 
pairs flow coherently in the presence of a scattering obstacle of arbitrary strength. The 
Josephson effect corresponds to the limit in which a strongly reflecting obstacle (typically, a 
tunneling barrier0 or a point contact!) reduces drastically the effective coupling between two 
bulk superconductors while still preserving global coherence. In the absence of an external 
bias, the current is given by the Josephson relation I = Jcsin(A^), where A(p is the phase 
difference between the two superconductors. The opposite limit is that of supercurrent flow 
in a perfect lead without appreciable scattering. In the appropiate gauge, this regime is 
characterized by a superconducting gap of uniform amplitude and a linearly varying phase 
whose gradient is proportional to the current. Specifically, in the Ginzburg-Landau limit, the 
current density can be written j = (eh/m)\tp\ 2 'Vip, where vp = \i\)\e lLp is the superconducting 
order parameter. 

An adequate measure of the scattering strength is the average transmission probability 
Tq for a Fermi electron passing through the barrier or contact in the normal phase, 



where is the device normal resistance, A is the cross section area of the semiinfinite leads, 
and kp is the Fermi wave vector. "Weak" and "strong" superconductivities are then charac- 
terized by T 1 and T ~ 1, respectively. For a structure in which superconductivity is not 
weakened by one-electron reflection, such as a S-N-S junction without current concentration, 
a more general parameter is Ic/Ib, where Ic is the critical current of the structure and Is 
is the bulk critical current of the perfect lead. It seems natural to ask how is the supercon- 
ducting flow for intermediate values of T Q or, more generally, Ic/Ib, i-e., how is the crossover 
between the two extreme limits of superconducting flow. This rather fundamental question 
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is of special current relevance, in view of the recent activity on superconducting mesoscopic 
structures (see, for instance, Refs.lhi). In the case of a superconducting point contact, the 
intermediate regime would correspond to contact widths not much smaller than the width 
of the semiinfmite leads. Alternatively, in the case of tunneling barriers, the crossover could 
be explored by considering different degrees of transparency at the Fermi level. In the case 
of a S-N-S junction, the intermediate behavior would be displayed by relatively thin normal 
metal layers located between two superconductors. 

A preliminary version of some of the results contained in this article has been briefly 
presented in Ref.i. 



II. SELF-CONSISTENCY AND CURRENT CONSERVATION 

Theoretical studies of weak superconductivity almost invariably assume that the phase 
is constant within the two superconductors. This is generally a reasonable approximation, 
since, by definition, in this regime, Ic "C Ib- As a consequence, the variation of the phase in 
the bulk of the superconductor displayed by current carrying solutions can be safely neglected 
in a wide range of length scales. It is clear that the approximation of an asymptotically 
uniform phase cannot be justified if Ic becomes comparable to Ib, which will certainly be 
the case in structures with moderate or negligibly weakened superconductivity. The more 
general situation will be that of a phase which varies linearly throughout the lead except in 
a finite region near the scattering center where it varies faster. 

In order to discuss some general questions related to self-consistency, we focus in this 
section on structures in which the decoupling between the two sides of a superconductor 
is due to one-electron scattering by a barrier or point contact. The conventional way of 
generalizing the BCS theory to the presence of an arbitrary one-electron potential is based 
on the Bogoliubov - de Gennes (BdG) equations@: 
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where H is the one-electron Hamiltonian, A is the gap function, and [u n (r) , v n (r)] and 
e n > are, respectively, the normalized wave function components and the energy of the 
quasiparticle n. The self-consistency condition for the gap function is3 

A = V£«X(l-2/ n ), (3) 

n 

where V is the electron-phonon coupling constant and /„ = [exp(e n /kT) + l] -1 . The BdG 
Hamiltonian can alternatively be written 

H = ~ Y C « dr K ; «( r )| 2 + Y tnlnvlno, (4) 
na na 

where 7^ CT creates quasiparticle n with spin a. In terms of the quasiparticle operators, the 
charge and current density operators are written 

P = e {Y Kf+Y (K U m ~ V* V m)lLlrna 
na nma 

+ U n Vm(T-fmaln,-a + U n V m°li,-allna} ( 5 ) 

nma nma 
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+ Y v m ~Du n a-f ma -f n ^ a - Y vlPKeTn-alLa}, ( 6 ) 

nma nma 

where e = — |e|, a = ±1, and D is defined as f~Dg = f(Vg)—(Vf)g. In Eqs. (5) and (6), the 
contributions from the condensate and the quasiparticles have been clearly separated. The 
quasiparticle contribution can in turn be divided into a part which conserves the quasiparticle 
number and a part which does not. The non- conserving components will not contribute to 
the expectation values < j > and < p > but will play an important role in the quantum 
fluctuations of the electronic charge and current densities. 

If one attempts to solve the BdG equations (2) in a given structure subject to the 
boundary condition that the phase takes certain constant values on each semiinfinite lead, 
one generally finds from (6) a nonzero value of the total current. This general feature can 
be illustrated by solving exactly a specific and very important example, namely, that of 
a strictly one-dimensional superconductor (i.e., with only one propagating channel for the 
Fermi electrons) with a barrier of arbitrary transmission Tq at the Fermi level. In this 



model, the phase is assumed to be uniform on each side of the barrier. A non self-consistent 
resolution of the BdG equations at zero temperature yields the current!!!'! 

T(/\ m \ - e l A l T sin (A^) 

{ V) 2h [l-T sin 2 (A^/2)]V2 ^ 

where Aip is the difference between the phases on each side of the barrier. Fig. 1 shows 
the current I (Aip) for several values of To. As the strength of the barrier decreases, the 
current departs from the ideal Josephson behavior and its maximum is displaced towards 7r. 
In particular, when T equals unity, the current is given by the formula!'! 

I(A<p) = (e\A\/h)sin(Aip/2), (8) 

with —7i < Aip < 7r and periodicity 2tt. This result is clearly not self-consistent, since a 
uniform phase should be associated with a vanishing equilibrium current, at least in the 
asymptotic region. Actually, a more detailed calculation revealsil that the current (6) is 
localized exponentially around the barrier in a region of width tt^/Tq sin(A</?/2), where 
£o = hvF/Tc\A\ is the zero-temperature coherence length. This peculiar feature can be 
traced back to the existence of a localized, current-carrying quasiparticle at the interfaced. 
Thus, one finds that the equilibrium current is nonzero near the scattering center and zero 
in the asymptotic region. In the steady state, this situation clearly involves a violation of 
charge conservation. Below we show that the relation between self-consistency and current 
conservation is in fact a general property of the BdG equations. 

The time derivative of the charge density operator can be computed by applying (4) and 
(5) to the relation p = (l/ih)[p, H\. The result is 
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+ {e n + e m )v[UnV m lmaln,-a + «7n -clLaW ( 9 ) 

which obviously yields < p >= 0, as expected from a stationary scattering description (we 
have used the properties < 7L7nV >= /n<W<W and < ^naln'a' >= 0). Combining (6) 
and (9) we obtain for the continuity equationll! 



<V-j> + <^>= x $: Im{A*« n <(l " 2/n)}. (10) 

By comparing this result with Eq. (3), it becomes clear that charge conservation is only 
guaranteed when the self-consistency condition is satisfied. In the language of Ref.0, the 
BCS-BdG theory is a conserving approximation only for solutions that satisfy the mean-field 
equations. It is interesting to note that the the condensate and quasiparticle contributions 
to the electric charge are not conserved separately, but only the sum of the two, and if 
the description is fully self-consistent. The relation between self-consistency and current 
conservation has also been noticed by Furusaki and TsukadaS, who have derived an equation 
similar to (10) in which condensate and quasiparticle contributions are however not clearly 
separated. This seems to lead to a misinterpretation. Unlike suggested in Ref.B, preservation 
of current conservation is not achieved in general by merely converting quasiparticle current 
into condensate currentSIll, but by truly implementing global self-consistency. A good proof 
of this assertion is that, within a non self-consistent scheme, the source term in Eq. (10) is 
generally nonzero even at zero temperature, when no quasiparticles exist. 

Before we proceed further a few additional remarks go in place. In one dimension, Eq. 
(7) is incorrect when T is not much smaller than unity. In particular, Eq. (8) is clearly 
wrong, since no bound quasiparticle should exist in the absence of a barrier. Of course, 
the main inconsistency lies in the very assumption of an existing phase difference, which 
cannot be maintained without a scattering obstacle (an abrupt change in the phase cannot 
survive the implementation of self-consistency). It will be seen in the following section that 
the appropiate generalization of the concept of phase difference to structures with arbitrary 
transparency is the phase offset, in terms of which the transparent limit will be quite different 
from (8). In studies of superconducting quantum point contacts, equations which generalize 
(7)i and (8)i'i to the presence of many transverse modes can be found. In these cases, 
the lack of formal self-consistency is justified. The localized nonzero current corresponds to 
the current in the vicinity of the point contact and the vanishing of the asymptotic current 
describes the widening of the contact into the reservoir. Therefore, Eqs. (7) and (8), as well 



as their multimode generalizationsHti, are correct as long as Ic Ib- This is the case when 
the number of propagating modes in the contact is much smaller than the number of modes 
in the wide leads. 



III. STUDY OF THE CROSSOVER 

From the discussion in the previous section, it is clear that, in order to achieve a uni- 
fied view of the crossover from weak to strong superconductivity, one must deal with self- 
consistent, current conserving solutions of the BdG equations in which a nonzero current is 
associated with a linearly varying asymptotic phase, and allow for arbitrary critical currents 
Ic — Ib- Unfortunately, the self-consistent resolution of the BdG equations for arbitrary 
currents is in general a demanding numerical task. By contrast, the formalism of Ginzburg- 
Landau (GL) provides a relatively simple method to learn about the global properties of 
those self-consistent solutions. Therefore, our goal in this section is to study the solutions 
of the GL equation for a one-dimensional superconductor in the presence of a delta poten- 
tial of arbitrary strength. Specifically, we wish to analyse the stationary solutions of the 
free-energy functional 

F = J dx[\V^\ 2 /K 2 - (1 - V 5(x)M 2 + H 4 /2] (11) 

where n = A/£ (A(T) is the penetration depth and £(T) is the temperature-dependent 
coherence length) and Abrikosov units are used. In these units, A(T) is the unit of length, 
the order parameter ip is measured in units of tpoo (absolute value of the bulk order parameter 
at zero current), and (ft/ra)( , 0^ o /£(T) is the unit for current. The complete crossover between 
weak and strong superconductivity will be explored by considering all values of the scattering 
strength g = kVq ranging from g very large (ideal Josephson behavior) to g = (uniform 
superconductor). In Eq. (11), F must be understood as the freen energy per unit area. This 
model should give a fairly adequate picture of a quasi-one- dimensional superconductor (of 
width w <C A, £) in which a (narrower) point contactB or a normal metal island has been 
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inserted. A clean point contact at low temperature could not be described by (11), since, 
in the weak superconductivity limit, this structure yields a current-phase relation of the 
type (8)IHIll instead of the usual sin(A<£>) behavior. On the other hand, the model (11) is 
not appropiate for a quantitative description of tunneling barriers because, in the limit of 
large g, the repulsive potential Vq5(x) yields hard-wall boundary conditions, which do not 
correspond to a GL description of the metal-insulator interface^. A similar model, with the 
5 function replaced by a square barrier, was studied by Jacobson0, who however focussed 
on the low current limit. Volkovll also used a delta function to describe a S-N-S junction 
but only analysed the small current case. 

If we factorize ip(x) = R{x)e %ip ^ x \ the GL equations take the form 

K- 2 d 2 R/dx 2 + [1 - V 8(x)](l - j 2 /R i )R - R 3 = 

dip/dx = nj/R 2 (12) 

where the current density j is a conserved number (I = jA). We are interested in solutions 
which satisfy the boundary conditions 

dR{x)/dx = 

(p(x) = qx ± A(p/2, for x — > ±oo (13) 

Current conservation requires the product R 2 (p' = nj to be constant, which can only be 
achieved with a nonzero q = nj/R 2 ^ in the asymptotic solution. The general solutions for 
R and if are of the form 

R 2 (x) = a + btanh 2 [Ku(xo + \x\)] 
<f(x) = qx + sgn(x){arctan[/9 tanh(K«(a;o + \x\))] 

— arctan[/9 tanh(KMX )]}. (14) 



In Eq. (9), o(2 - a) 2 = 8j 2 , with < a < 2/3, 6=1- 3a/2, u = yJb/2, (3 = ^Jb/a, and x 
is obtained from the matching condition at the site of the delta potential, which gives rise 
to the cubic equation 
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a/26/3i/o(1 - vl) - 9(1 + PVo) = 0, (15) 

where yo = tanh(Kuxo) and thus only the solutions satisfying < yo < 1 are of interest. 
The solutions turn out to be uniquely parametrized by the phase offset Aip, whose general 
expression is 

Aip = 2[arctan(/3) — arctan(/3yo)]- (16) 

The resulting curve j(Aip) is displayed in Fig. 2. The inset shows the critical current 
as a function of g. The Josephson limit is well achieved for g > 8 while jc saturates 
to j B = 2/3\/3 = 0.385 as g — > 0. For large g, one finds the ideal Josephson behavior, 
j = jcsin(Aip), with jc = l/2g taking small values. For g = 0, two types of solutions are 
obtained. One of them is entirely expected: for Aip = 0, all currents are possible ranging 
from j = to j = 2b- These are the solutions of the uniform superconductor in which ip' 
and R take constant values. The second type of solutions are the solitons of the ip A theory 
defined by (11) for arbitrary values of the current ja. These kinks separate two domains in 
which the phase varies linearly, 

<f(x) = qx + arctan[/3 teuah^nux)} (17) 

with a total phase offset of Aip = 2arctan(/3). It is interesting to note that, unlike in the 
j = case, the phase offset (which here plays the role of the soliton charge) can be different 
from 7r. These solitonic solutions are equivalent to the saddle-point configurations which 
were considered by Langer and Ambegaokar0 in their study of the resistive behavior of 
one-dimensional superconductors. 

Fig. 3 shows the behavior of R(x) and <p{x) for g = 0.2 and two different values of the 
current. In Fig. 3a, it is clearly seen that, for j = 0.01, the solitonic solution almost vanishes 
at x — 0. For the same solution, Fig. 3b shows that the spatial variation of (p(x) is almost 
negligible except for a step-like feature at x = (the phase can be shown to vary in a length 
scale j/n if j is small). For j = 0.35 (close to js), the phase displays a linear increase with 
x with an offset due to a faster variation in the vicinity of x = 0. 

9 



An interesting feature of the j(Aip) curves which can be clearly observed in Fig. 2 is that, 
as the scattering is turned off, the maximum current is displaced towards lower values of 
Aip. This is in sharp contrast with the behavior shown in Fig. 1 for the non self-consistent 
solutions. It has already been noticed that a superconducting point contact displays the 
same behavior as its propagating channels evolve from low to high transmission^. There is 
of course no contradiction between our results and those obtained for point contacts, since 
the latter apply only in the limit Ic <C Ib, while the low g curves in Fig. 2 are only relevant 
in the Ic ~ Ib case. 

Baratoff et a/.0 considered a S-S'-S structure in which S and S' are two dirty super- 
conductors of differing properties. As a function of the similarity between S and 5", they 
obtained results which qualitatively resemble those obtained by us. However, their focus 
was not in the crossover from weak to strong superconductivity, but rather in the quali- 
tative modeling of weak links. In particular, they did not consider the S = S' case and, 
although the relation to Ref.0 is noticed, no association is made between the branch to the 
left of the maximum in the j—A(p curve and the trivial solutions of the uniform case. More 
recently, Kupriyanovlil has studied the properties of an S-I-S structure by means of the 
Usadel equations, which apply in the dirty limit. He considers several values of the barrier 
transparency and obtains results which, after a nontrivial scale transformation (the phase 
change across the barrier instead of the phase offset is used as a parameter) , can be shown 
to be qualitatively similar to those displayed in Fig. (2). However, in the transparent limit, 
no mention is made in Ref.0 of the relation to the solitonic solutions of Ref.0 nor to the 
uniform solutions, as discussed here by us. 

IV. CROSSOVER IN LONG BRIDGES. BREAKDOWN OF THE AC 

JOSEPHSON EFFECT. 

So far we have focussed on the relation between the current j and the phase offset Aip, 
which uniquely parametrizes the solutions of the GL equations (12). However, it is also 
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convenient to plot the current as a function of the total phase difference x between two 
reference points. These two points can be, for instance, the extremes of a superconductor 
of length L with an effective barrier in its center. A typical case would be that of a narrow 
bridge connecting two wide reservoirs through smooth contacts beyond which the phase 
gradient can be safely neglected!. For a given length L, one can compute x horn the 
relation x — l p(L/2) — ip(—L/2). If L ^> (ku)~ 1 , we can approximate 

X^qL + Aif. (18) 

Since u — > as j — > ja, there is a threshold current jth(L) above which Eq. (18) does not 
apply. For kL ^> 1, jth — 3b(1 — 27/16k 2 L 2 ). In Fig. 4, the resulting curve j(x) is plotted 
for kL = 10. In such a case (Jb ~ 3th) I 3b — 0.016. It can be observed that, for large g, the 
ideal Josephson behavior is displayed, while, for sufficiently small g, the current becomes a 
multivalued function of the total phase x- The pattern shown in Fig. 4 is actually repeated 
periodically with a period of 2tt. In the case of g small it becomes clear from the comparison 
with Fig. 2 that the upper branch corresponds to solutions with a linearly varying phase 
(Aip ~ 0), while the lower branch is given by the solitonic solutions with a nonzero phase 
offset. This feature has also been noticed recently by Martin- Rodero et oi@, who have 
computed numerically the self-consistent solutions of the BdG equations for a linear chain 
coupled to two Bethe lattices at zero temperature. The discontinuity in the derivative at 
the top of the g = curve in Fig. 4 reflects the discontinuous transition from the uniform 
to the solitonic branch shown in Fig. 2. However, this cusp cannot be observed in bridges 
of finite length since it always lies above the threshold of validity of Eq. (18). 

In Fig. 5 we display the phase of the order parameter as a function of the position and the 
phase offset for g = 10 (Josephson limit). When A(p = ±7r, the current is zero. This requires 
an abrupt jump of ±7r at x — 0, which is possible because R(0) = in these solutions, as 
can be proven quite generally. These are the phase-slip configurations which permit the 
existence of the ac Josephson effect. As an external driving voltage is applied between two 
points on different sides of the junction, the phase is forced to vary at a constant rate and 
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the whole sytem responds adiabatically by evolving along the continuous set of stationary 
solutions. The existence of these step function solutions makes it topologically possible for 
the phase at every point to increase both monotonically and continuously with time. Since 
R(0) = 0, the two superconductors are completely decoupled and Aip = n is equivalent 
to Aip = —7i. As the system is driven by the external bias through the different values 
of Aip and reaches the value Aip = it, it automatically reenters through the topologically 
equivalent configuration with Aip = —n and the phase at the boundary can continue to 
increase monotonically. Thus the existence of the ac Josephson effect relies on the ability 
of the system to undergo adiabatic phase-slips under the action of an external bias. It is 
interesting to note that, at the particular value Aip = ±7r, the configuration of the order 
parameter is independent of g, since then i?(0) = 0. In particular, it is identical to the 
phase-slip configuration in the absence of a barrier, as studied by other authors (see, for 
example, Refsffl). 

At zero current, all points on one side of the barrier have the same phase. In particular, 
ip(x > 0) = ±7r/2 for Aip = ±7T. By contrast, the solutions with nonzero current have an 
asymptotic phase which grows linearly with position, as shown in Eq. (17). Thus, for suffi- 
ciently large x, it is not possible to have (p(x) increasing monotonically as Aip varies between 
the two equivalent configurations with Aip = ±7r. As a consequence, the system cannot re- 
spond adiabatically to a constant voltage being applied at points that are sufficiently distant 
from the junction. The only choice for the system will be to undergo nonadiabatic, fluctu- 
ating processes of the type studied by Langer and AmbegaokarS (albeit with g 0), which 
will originate a resistive behavior. The threshold for this type of response is given by the 
condition 

dip(x b ) 



= at Ay? = 7T (19) 

Al/?=7T 



dAip 

If one identifies Xj, = L/2 and <p(xb) = x/2, this is also the condition for the onset of 
bivaluedness in the j(x) curve of Fig. 4, which requires dx/dj = at x — ^ (note that, 
dj /dAip 7^ at Aip = n). Thus, if the electrodes are applied at points \x\ > Xj,, there is 
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a breakdown of the ac Josephson effect due to the fundamental inability of the system to 
respond adiabatically to that particular type of external constraint. 

Let us estimate the breakdown length x^. For g large, one can show that Xb — g/^- We 
notice at this point that the value of the parameter g can be adjusted to a realistic setup 
by exploiting the relation 

g = 1.30 -(jb/jc), (20) 

which applies in the Josephson regime, and noting that 3b 1 3c = Ib/Ic- We have considered 
explicitly four types of structures which are known to display a standard sin(Ayj) behavior 
in the Josephson limit for T close to T c : (a) a tunnel junction with average transmission T 
for the Fermi electrons, (b) a clean point contact with average transmission T , (c) a narrow 
bridge between two superconductors made of a dirty normal metal of length L and coherence 
length £jv at T ~ T c , and (d) a S-N-S structure without current concentration (N and S 
have the same width). Cases (a) and (b) fall within the same category in the GL limit, with 
an expression Iq = n /S? {T) / AeR^ksT for the critical currentll. Noting that, for T close to 
T c , the gap function and the order parameter are related by0 ip = 0.326 y/nA/kBT c , where 
n is the electron number density, we arrive at 

g- 1 ~2.0T (Z(T)/Z ). (21) 

For case (c), the critical current is! I c = {AA 2 (T) /ixeR N k B T)(L/ £, N ) exp(-L/£ N ), if L > 
£tv. Thus one obtains 

g- 1 ~ 3.23 To(£(T)/£o)(L/&v) exp(-L/^). (22) 

For a S-N-S structure without current concentration, the critical current isH Ic = 
A(ehn/2m)(\T — T C \/T C )(^ N /^ 2 (T)) exp(— L/^ N ). As a consequence, 

g- 1 ~ 10.6 (W£C0) exp(-V6v). (23) 

Shifting to real units, we arrive at the relations 

13 



x b ~ 0.50 to/T 



(a),(b) 



x b ~ 0.31 (^ N /T L)e L/ ^ 
x b ~ 0.094 (£ 2 (T)/£v)e L /« 



(d) 



(c) 



(24) 



for the maximum distance at which a constant voltage can be applied in order to observe 
the ac Josephson effect. 



We have studied the nature of the crossover from ideal Josephson behavior between two 
weakly coupled superconductors to bulk superconducting flow in a perfect superconducting 
lead. We have argued that a self-consistent resolution of the BdG equations is mandatory 
in a microscopic study of the crossover and have proved that charge conservation is only 
guaranteed when the requirement of self-consistency is satisfied. We have performed a study 
of the crossover by solving exactly the Ginzburg-Landau equation for a one-dimensional 
superconductor in the presence of a delta potential of arbitrary strength. The pairs of 
Josephson solutions with equal current have their scattering free counterparts in the pairs 
formed by a uniform and a solitonic solution. This relation has allowed us to understand 
some aspects of the multivalued current-phase relation in narrow bridges. The complete 
knowledge of the set of stationary solutions for different values of the scattering strength g 
has helped us to gain a more detailed understanding of the adiabatic response to a constant 
external bias, which has been shown to rely on the feasibility of adiabatic phase-slips. If a 
voltage is applied at points which are sufficiently far from the junction, there is a breakdown 
of the Josephson effect due to the intrinsic impossibility of changing adiabatically the phase 
at a distant point in a continuous and monotonic manner. 



V. CONCLUSIONS 
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FIGURES 

FIG. 1. Current j as function of the phase offset A 92 for a non-self-consistent solution The 
curves labeled a, b, c, and d are for the cases To = 0.999, 0.99, 0.9 and 0.4, respectively. 

FIG. 2. Current j as function of the phase offset A93. The curves are labeled a, b, c and d for 
the cases g = 0, 0.5, 3 and 10, respectively. Inset: critical current jc versus scattering strength g; 
solid line gives the the exact result and dotted line corresponds to the Josephson limit l/2g. 

FIG. 3. The amplitude (a) and the phase (b) of the order parameter plotted as a function of 
position in the g = 0.2 case (re = 1), for values of the current j = 0.01 (curves labeled a) and 
j = 0.35 (curves labeled 6). 

FIG. 4. Same as Fig. 2, for the total phase difference \ between the extremes of a supercon- 
ductor of length L = 10. 

FIG. 5. The phase of the order parameter is plotted as a function of the position x (re = 1) and 
the phase offset Aip for the case g = 10. 



18 



